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Abstract 

We investigate the large-time behavior of solutions to an outflow problem of the full compressible 
Navier-Stokes equations in the half line. The non-degenerate stationary solution is shown to be 
asymptotically stable under large initial perturbation with no restriction on the adiabatic exponent 
7, provided that the boundary strength is sufficiently small. The proofs are based on the standard 
energy method and the crucial step is to obtain positive lower and upper bounds of the density and 
the temperature uniformly in time and space. 
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1 Introduction 

We study an initial-boundary value problem for the full compressible Navier-Stokes equations 

! Pt + ( pu) x = 0, 

(pu) t +(pu 2 + P) x = (pu x ) x , (1.1) 

(, pE) t + (puE + uP) x = ( k9 x + puu x ) x , 

where t > 0 is the time variable, x G R+ := (0,oo) is the spatial variable, and the primary dependent 
variables are the density p , fluid velocity u and absolute temperature 0. The specific total energy E = 
e + ^ it 2 with e being the specific internal energy. The viscosity coefficient p and the heat conductivity k 
are assumed to be positive constants. Here we focus on the ideal polytropic gas, that is, the pressure P 
and the specific internal energy e are given by the constitutive relations 

P = Rpd , e = c v 0, (1.2) 

where R > 0 is the gas constant and the specific heat c v = R/{ 7 — 1) with 7 > 1 being the adiabatic 
exponent. 

The system (1.1) is supplemented with the initial data 

(p, u, 0)|t=o = (Po, u 0 , do), (1-3) 
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which are assumed to satisfy 

lim (p 0 ,u 0 ,9 0 )(x) = (p+,u+,9 + ), (1.4) 

x—too 

where p+ > 0, u+ and 0+ > 0 are constants. For boundary conditions, we take 

u(i,0) = u_, 9(t, 0) = 9_ > 0, (1.5) 


where u- < 0 and 0 _ are constants. 

The assumption u— < 0 means that the fluid blows out from the boundary, and hence the problem 
(1.1)-(1.5) is called the outflow problem. In the case of u- = 0 the problem is called the impermeable 
wall problem, while the problem (1.1)-(1.5) with the additional boundary condition p(t, 0) = p_ in the 
case of U- > 0 is called the inflow problem. 

Matsumura [15] considered initial-boundary value problems for the isentropic Navier-Stokes equations 
in the half line R+ and proposed a complete classification about the precise description of the large¬ 
time behaviors of solutions. Since then, some results have been obtained for the rigorous mathematical 
justification of this classification. For the impermeable wall problem, the stability for the viscous shock 
wave is obtained in [16] for small initial perturbation, while Matsumura-Nishihara [17] showed that the 
rarefaction wave is asymptotically stable under large initial perturbation. For the inflow problem under 
small initial perturbation, see [18] for the stability of the boundary layer solution and its superposition 
with the rarefaction wave, and see [3] for the stability of the viscous shock wave and its superposition 
with the boundary layer solution. Recently, Fan et al. [1] established the asymptotic stability of both the 
boundary layer solution and the supersonic rarefaction wave for a certain class of large initial perturbation. 
For the outflow problem, the stability of the stationary solution and its superposition with the rarefaction 
wave is proved in [9, 10] under small initial perturbation, while a convergence rate of solutions toward the 
stationary solution was obtained in [19] provided that the initial perturbation belongs to some weighted 
Sobolev space. Huang-Qin [5] improved the results in [9, 10] to large initial perturbation. 

There have been some works on the large-time behaviour of solutions to the initial-boundary value 
problems of the full compressible Navier-Stokes equations (1.1)-(1.2) in the half line. For the inflow 
problem, Qin-Wang [23, 24] showed the existence of the boundary layer solution as well as the stability 
of the boundary layer and its superposition with the viscous contact wave and rarefaction waves. The 
asymptotic stability of the rarefaction wave, boundary layer solution, and their superposition is proved 
in [2] for both the impermeable wall problem and the inflow problem. Kawashima et al. [ 8 ] proved 
the existence, the nonlinear stability and the convergence rate of the stationary solution for the outflow 
problem. We note that the results in [2, 8 , 23, 24] are all concerned with small initial perturbation. 
Thus a problem of interest is whether stability results on the stationary solution and the rarefaction 
wave hold for the outflow problem (1.1)-(1.5) with large initial perturbation as in [5] for the isentropic 
model. In this direction, motivated by [21], Qin [22] proved that the non-degenerate stationary solution 
is asymptotically stable under “partially” large initial perturbation with the technical condition that the 
adiabatic exponent 7 is close to 1 . 

In this article, we establish the large-time behavior of solutions toward stationary solutions to the 
outflow problem (1.1)-(1.5) under large initial perturbation without any restriction on the adiabatic 
exponent 7 . We expect that the solution converges to a stationary solution ( p,u,9)(x ) of (1.1): 

(pu)' = 0, ' := zeR+, 

| (pu 2 + Py = pu ", (i-®) 

[ (puE + uP)' = u9" + p(uu ')', 

where P := Rp9 and E := c v 9 + \u 2 . We also assume that the stationary solution satisfies the far-field 
condition (1.4) and the boundary condition (1.5): 

(u, 9)( 0) = (u-, 0_), lim (p, u, 9)(x) = ( p+,u+,9 + ). (1.7) 

x—too 

We define the boundary strength 6 as 


S := |(u+ - U-,9+ - 6 L)|, 
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and the Mach number at infinity as 


M+ 


l u +l 

* 

c+ 


where c+ := \JR"/9 + is the sound speed. The existence and the properties of the stationary solution 
(p,u,0) satisfying (1.6) and (1.7) are quoted in the following lemma. 


Lemma 1.1 (Existence of stationary solution [8]). Suppose that ( U-,9 _) satisfies 


(u_, 0_) G M + := {(it, 9) G R 2 : |(u — u+, 9 - 6 + )\ < <5 0 } 


( 1 . 8 ) 


for a certain positive constant <5o. 

(i) For the case M + > 1, there exists a unique smooth solution ( p,u,9) to the problem (1.6)-(1.7) 
satisfying 

\d x {p~ P+,u- u+,9 - 9+){x)\ < CSe~ cx (1.9) 

for integer n = 0,1, 2, • • ■ , where C and c are positive constants. 

(ii) For the case M + = 1, there exists a certain region C such that if (it_,0_) G Ai°, then 

there exists a unique smooth solution ( p,u,6) to (1.6)-(1.7) satisfying 

rjxn+i 

I TO - P+, u u + ,9 - 9 + )(x) | < + +1 + C6e~ cx (1.10) 

for integer n = 0 , 1 , 2 , • • • . 

(iii) For the case M + < 1, there exists a curve AA~ C AA + such that if(u-,9_) G Al + , then there exists 
a unique smooth solution ( p,u,9 ) to the problem (1.6)-(1.7) satisfying (1.9). 

Our main result is now stated as follows. 

Theorem 1 . Assume that there exists an stationary solution (p,u,9) to the problem (1.6)-(1.7) satisfying 
(1.9). If 

inf (po, 9o)(x) > 0, (p 0 - p,u 0 - u,9 0 - 9) G H 1 (R+), (1.11) 

then there is a positive constant e\ such that if S < ei, then the outflow problem (1.1)-(1.5) admits a 

unique solution ( p , it, 9) satisfying 


(p — p,u — u,9 — 9) G C([0,oo); H 1 ( R+)), 

Px-Px G L 2 (0, oo;L 2 (K+)), (u x - u x , 9 X - 9 X ) G L 2 ( 0, oo; Jf 1 (R + )), 


( 1 . 12 ) 


and 

lim sup \{p— p, u — u, 9— 9)(t,x)\ = 0. (1.13) 

t - K ’ 0 xeR+ 

Remark 1.1. Theorem 1 shows that the non-degenerate stationary solution (i.e. the one which decays 
exponentially as (1.9)) is asymptotically stable even for large initial perturbation with general adiabatic 
exponent 7 , provided that the boundary strength is small. 

To derive the large-time behavior of solutions toward stationary solutions, it is sufficient to deduce 
certain uniform (with respect to the time t) a priori estimates on the perturbations. The essential step 
is to get the lower and upper positive bounds for both the density p and the temperature 9 uniformly in 
time and space. In case of small initial perturbation, as in [8], we can use the smallness of the a priori 
H 1 -norm of the perturbation to obtain the uniform bounds of the density p and the temperature 9. Based 
on such uniform bounds and the smallness of the boundary strength 6, one can derive certain uniform a 
priori ii/^-norm energy type estimates on the perturbation. In case that the adiabatic exponent 7 is close 
to 1, as in [22], we can control the lower and upper bounds of 9, and by using the Kanel’s technique, we 
can then obtain the uniform positive bounds for p provided the boundary strength 8 is small. 

We are interested in showing the stability of the non-degenerate stationary solution under large initial 
perturbation with general 7 . For this purpose, we first deduce the basic energy estimate with the aid 
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of the decay property of the non-degenerate stationary solution provided that the boundary strength S 
multiplied with a certain function of the a priori lower and upper bounds of density p and temperature 9 
is suitably small (see Lemma 2.2). Next, to get uniform pointwise bounds of the density p, we transform 
the outflow problem (1.1)-(1.5) into a free boundary problem in the Lagrangian coordinate. By modifying 
Jiang’s argument in [ 6 , 7] for fixed domains, we will use a cut-off function with parameter to localize the 
free boundary problem, and then we will deduce a local representation of the specific volume v = 1/p 
to establish the uniform bounds of v. With such uniform bounds of the density p in hand, we can 
derive the iL 1 -norm (in the spatial variable x) estimate of the perturbation 9(t,x ) — 9(x) uniformly in 
the time t in the Eulerian coordinate. Then the uniform upper bound of the temperature 9 is obtained 
in light of the Cauchy’s inequality. We note that our derivation of the uniform upper bound of 9 is 
motivated by the recent work [14], which is concerned with the stability of the constant state to the 
compressible Navier-Stokes equations on unbounded but fixed domains in the Lagrangian coordinate. 
The maximum principle gives us the lower bound of the temperature 9 locally in time t. In view of the 
a priori assumption ( 2 . 6 ), we have to get the uniform positive lower bound of the temperature 6 , which 
will be achieved by combining the local lower bound of 9 and the detailed continuation argument. 

Another interesting problem is on the asymptotical stability of the rarefaction wave and its superpo¬ 
sition with the non-degenerate stationary solution to the outflow problem (1.1)-(1.5) under large initial 
perturbation. See [5, 10, 11] for the isentropic case and [22] for the case of small initial perturbation. We 
will consider this problem in a forthcoming paper. 

The layout of this paper is as follows. After stating the notations, in section 2, we establish the 
desired energy estimates and obtain the uniform bounds of both density and temperature. In section 3, 
we extend the local solution step by step to a global one and prove the stability of the stationary solution 
by combining the a priori estimates obtained in section 2 with the the continuation argument. 

Notation. We employ C or Ci ( i £ N) to denote a generic positive constant which is independent of t, 
x and S. Notice that all the constants may change from line to line. The Gaussian bracket [x] means 
the largest integer not greater than x. For function spaces, L P (M. + ) (1 < p < oo) stands for the usual 
Lebesgue space on R + equipped with the norm ||-||lp and H k ( R + ) the usual Sobolev space in the L 2 
sense with norm || • ||fc. We use the notation || • || = || • || L 2 . We denote by C k {I]H p ) the space of fc-times 
continuously differentiable functions on the interval I with values in H p { R + ) and L 2 (I;H P ) the space 
of L 2 -functions on / with values in H p ( R+). 

2 A priori estimates 

2.1 Reformation of the problem 

This section is devoted to deriving a priori estimates on the solution (p, u, 9) to the outflow problem 
(1.1)-(1.5). To this end, we regard the solution ( p,u,9 ) as a perturbation from the stationary solution 
(p, u, 9) and put the perturbation )(p, ip, $) by 

(0, ip, $)(t, x) := (p, u, 9)(t, x) - (p, u , 9)(x). (2.1) 

Subtracting (1.6)-(1.7) from (1.1)-(1.5) yields 

(,At + u(p x + pip x = -u x (p - p x ip, t > 0, x £ K + , 

p(ipt + uip x ) + (P - P) x = pip xx + p, 

< c v p('d t + u'd x ) + P4> x = kii) xx + pi/l +h, ( 2 . 2 ) 

(</>,-fM)|i=o = 

(tM)U=o = ( 0 , 0 ), 

g := -u x {u(f) + pi/), h := -c v 9 x (u(j) + pip) - u x (P - P) + 2pu x ip x , 


where 


(2.3) 




Outflow problem for full compressible NS 


5 


and the initial condition (cp 0 , ipo,&o) := (Po — p, ug — u, 9 0 — 9) satisfies 

((po, ipo, d 0 )(x) -4 (0,0,0), as a; —>■ oo. (2.4) 

The solution space X(0, T; mi, Mi; m 2 , M 2 ) is defined by 

X(0,T;mi,M i; m 2 ,M 2 ) := {(<p,ip,0) £ C([0,T]; M 1 ) : £ T 2 (0.Tji/ 1 ), 

(p x £ T 2 (0, T; T 2 ), mi < <j> + p < Mi, m 2 < i? + 9 < M 2 }. 

We summarize the local existence of solutions to the problem (2.2) in the following proposition, which 
can be proved by the standard iteration method (see [4]). 

Proposition 2.1 (Local existence). Suppose that the conditions in Theorem 1 hold. If ||(^o> ipo, 'do)||i < 

M, Ai < (f>o(x) + p(x) < Ai, and A 2 < dg(x) +9(x) < A 2 hold for all x £ R+, then (2.2) admits a unique 
solution (cp,ip,d) £ X (0, To; |Ai, 2Ai; |A 2 ,2A 2 ) for some constant T 0 = T 0 (Ai,A 2 ,M) > 0 depending 
only on Ai, A 2 and M. 

We now turn to deduce certain uniform a priori estimates for the perturbation (<p, ip, 9) £ X(0,T; mi, M\, m 2 , M 2 ) 
in the Sobolev space H 1 . Before doing so, we recall that C or C, (i £ N) will be used to denote some 
generic positive constant which depends only on mfR + (p 0 ,do) and ||(</>o>' ( ('0)'do)||i- For notational simplic¬ 
ity, we introduce A < B if A <CB holds uniformly for some constant C depending only on inf R+ (p Q , 9g) 
and ||(</. 0 ,^o,do)||i. The notation A ~ B means that both A< B and B < A. Besides, we will use the 
notation (p, 9) = (</> + p, d + 9) so that 

mi < p(t,x) < Mi, m 2 < 9{t,x) < M 2 , for all (t, x) £ [0,T] x R + . (2-5) 

Without loss of generality, we may assume that m, < 1 < Mi for i = 1,2. 


2.2 Basic energy estimate 

First, we have the following basic energy estimate. 

Lemma 2.2. There exists a sufficiently small eg > 0 such that if 


l(mi, Mi, m 2 , M 2 )<5 < eo, a(mi, Mi, m 2 , M 2 ) := m 1 u M 1 u m 2 U M 2 , 


then 


where 


/ + + / 4(«.°)*+/ / + 

J R + P Jo P Jo JR + P 

[ ps+ r P <s>fpy s ,o)ds+ f [ 

J r + Jo \PJ Jo Jr. 


€ , C < 

9 9 2 ’ 


A , 1 


£ := R9$ - + ^ +Cv 9<f> = , $(z) :=z-lnz-l. 


,PJ 2 

Proof. By a straightforward calculation, we find 


(p£)t + 


pu£ + ip(P - P) - /.iipifx - k-A- 

U 


9ip 2 x 9d 2 x 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 


1Z = p9 x 


Ru<f> 


+ c v u& ( — — c v ip In 

1 6 


— Rip In 


- pu x ip 2 


^ C v puO x + P'Ux'j 


(<pd d 2 \ u x ip x d 9 x dd x u xx <pip 

I —- + -=—14- All. -1- K,-t- ll - 


\p9 99 


2/i 


+ K~ 


9 2 


where 
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The identities 


and 


$(z)= f f e 1 ^\i + e 1 e 2 ( z -i))de 2 de 1 {z-i) 2 
Jo Jo 

1-1 (z - l)d0! 


In 0 = 


/ 


1 + 0\ {z — 1 ) 


imply 


$(*) + | In z\ 2 < (z- 1 + 1 ) 2 (z - l) 2 , (z + 1 )~ 2 (z - l) 2 < $(*). 
In view of (2.5) and (2.11), we apply Cauchy’s inequality to 1Z and obtain 


( 2 . 11 ) 


n< 


(px, &X1 8 X , 


AfimiV + Mito )" 2 # 2 + M i^> 2 + 1 f + 


2 , 

6 > 6> 2 


Integrate (2.10) over [0, T] x K+ and use the good sign of u- to derive 


p£ + ( - ) (s, 0 )ds + 


/R+ 


0 J R+ 


< i + 


/ j [px^'d x ,6 x , : u xx j 

Jo JR, V 7 


8 8 2 

.- 2/2 


Mim^ z (j) z + Mimj 2 )? 2 + M^ 2 + ^ + 


1.2 , V>* , ^ 


0 6» 2 


( 2 . 12 ) 


To estimate the terms on the right-hand side of (2.12), we use an idea in Nikkuni-Kawashima [20], i.e. a 
Poincare type inequality 

|</?(M)I < MM)| + y/x\\ip x (t)\\, x £ R+. 

Applying this inequality to <j>, if) and ■&, we deduce from (1.9) and (2.11) that 


J2 f [ d x(P,u,6) \d l x ((/),ip,d)\' 

_n i J 0 J R-l 


1 = 0,1 


< 6 


/o L 


i < / > ( s >o)i 2 + \\(^x,i>x,^x) (s)ir 


< M 2 m 1 1 S yo<J> [ - ) (s, 0)ds 
Jo \PJ 

+ 8 [ ||0x(s)|| 2 ds + M 2 8 f ( 

Jo Jo Jr 


0 JR + L 


ds 


e e 2 


(2.13) 


for each integer k = 1, 2, • • • . Plug (2.13) into (2.12) to derive 


p£+ p $ 


(s, 0)ds+ f f 
Jo J Rj 


6 6 2 


< 1 + m x i Mfm 2 2 8 j ( - ) (s, 0)ds 


+ m 1 2 Mim 2 2 8 f \\(f> x (s)\\ 2 ds + m 1 2 Mim 2 2 M 2 f f 

Jo ~ Jo Jr 


which implies 


f p£ + f p$ f-^J (s, 0)ds + f f 
J r + Jo \Pj Jo Jr. 

<l+m^ 2 Mim 2 2 8 f ||^ x (s)|| 2 (is, 
Jo 


0 JR + L 


^ , € 
e e 2 


e e 2 


(2.14) 
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provided (2.6) holds for a suitable small eo > 0. 

To control the last term of (2.14), we differentiate (2.2) 1 (first equation of (2.2)) with respect to x, 
and then multiply the resulting identity and (2.2) 2 by % and respectively, to discover 


' <t>l . 


'u<t>l 


Pxip 2 

. RHl 

.2 p 3 pp . 

t 

2 p 3 

pp 

pp _ 

x PP 2 


with 


Q = i%_ 

P HP 


2 Pxi> 

pp pp. 


uu x R9 X 

r> ~ 1 " Q 

PP 

Ux 


PP 
Ipx - 


^XX 

,3 


P 

_ PUx 
pp pp 2 


= Q 


2 Px^x . RPx 9 


(2.15) 


Pxx 

n 3 


- 


p 0 pp z 

p x u x (t>4> Pxx4> 2 


pp 


pp 


We integrate (2.15) over [0, t] x R + and use Cauchy’s inequality to have 

p^+ rr \qi 

U Jo lli 


7 +7 » 4 <x 

S + P Jo P Jo J R + P 


(2.16) 


Cauchy’s inequality yields the bound 


|Q|<^ + C(e) 


i9 2 


+ | (px i Ux, 9 X j Pxx: Hxx') 11 (<7 9 1 4*x j Ipx ) | • 


e p 2 

Plugging this inequality into (2.16), we take e > 0 suitable small and utilize (2.13) and (2.14) to obtain 


/R. P' 


P‘ 


/% s .o )+ /7 

Jo P Jo j R-; 

$2-1 n ft 

ip 2 + ~K +m 1 3 Mim 2 2 M%5 / ||<() a; (s)|| 2 ds 

Jo 


< 1 + 


/ 0 J M_|_ 

4 71 /r2 71/T2 


'0 

C 

9 

/ (p) ^ s ’° S > ds + J 


H ^ 

0 (9 2 


< 1 + 


7* 


tf 2 

X 

6 > 


- m7 10 M}°mo W Mi 0 5 


'0 «/]&+ L 

Taking eo > 0 small enough, we can derive from (2.6) that 

a2 N N r OA 2 


/ f + />»> + /7 

J R + P Jo P Jo J R+ P Jo J Rj 


Q2 "I 


7x 


(2.17) 


The estimate (2.7) follows by plugging (2.14) into (2.17) under the condition (2.6) for a sufficiently small 
eo > 0. We plug (2.7) into (2.14) to deduce the estimate (2.8). The proof of the lemma is completed. □ 


2.3 Pointwise estimates of density 

Once the basic energy estimate (2.8) is obtained, we can proceed to deduce the positive lower and 
upper bounds of the density uniformly in time. For this purpose, it is easier to consider the problem in 
the Lagrangian coordinate than in the Eulerian coordinate. We introduce the Lagrangian variable 

p(s,0)ds+ f p(t,z)dz, (2-18) 

Jo 

and define (p, u, 9)(t , y) := (p, u, 9)(t, x ). By the coordinate change (t, x) i-)- (t, y), the domain [0, T] x R + 
is mapped into 

:= {(t,y) : 0 <t<T,y> Y(t)}, 
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where 


Y{t):=—u- I p(s,0)ds. 
Jo 


The outflow problem (1.1)-(1.5) can be transformed into the problem in the Lagrangian coordinate: 


v t -u y = 0, 

Ut + Py = (^A 




CyO ~\~ 


+ ( P u )v — 


V Sy 

K,9y /MlUy\ 

V V ) 


V 


CM) \y=Y(t) = 

(v,u, 9)\t=o = (v 0 ,u 0 ,9 0 ), 


V > Y(t), 


(2.19) 


where v = 1/p stands for the specific volume of the gas. We drop the hats in the formula in this subsection 
for simplicity of notation. The basic energy estimate (2.8) in Eulerian coordinate can be transformed 
into a corresponding estimate in Lagrangian coordinate as a corollary of Lemma 2.2. 

Corollary 2.3. If (2.6) holds for a sufficiently small eg > 0, then 


sup 

o <t<T 



d 2 

u v 

v9 2 


< 1 . 


( 2 . 20 ) 


Note that the function Y(t) describing the boundary in the Lagrangian coordinate is part of the 
unknown, i.e. the problem (2.19) is a free boundary problem. To obtain the uniform bounds of the 
specific volume v for the free boundary problem (2.19), we introduce 


Qi(t) 


[Y(t),[Y(t)} + 2], i=[Y(t)] + 1, 

[i,i + 1], i > [y( t )] + 2 


( 2 . 21 ) 


for any integer i > [F(f)] + 1. Based on the basic energy estimate (2.20), we have the following lemma. 

Lemma 2.4. Suppose that (2.6) holds for a sufficiently small eg > 0. Then there exists a constant Cg > 0, 
which depends only on infn + (po,$o) and ||(<(>o, V'Oi'$o)||i; such that for all pair ( s,t ) with 0 < s < t < T 
and integer i > [Y (f)] + 1, 


C 0 1 < f v(s,y)dy, f 9{s,y)dy <Cg, (2.22) 

and there are points ai(s,t),bi(s,t) £ f h(t) such that 

Cg 1 < v(s,ai(s,t)),9(s,bi(s,t)) < C 0 . (2.23) 

Proof. Let 0 < s < t < T and i > [E(t)] + 1. From the definition of Y(t) and the sign of u we have 
Y(s) < Y(t ) and C [F(s),oo). In view of (2.20), we get 

[ $ (v) (s,y)dy+ ( $ (V) (s,y)dy < 1 . 

Jn.it) Jn.it) \9j 

Apply Jensen’s inequality to the convex function $ to obtain 

<3> f ToTTVT / s iV) d y ) ( lomi / §( s ^2/)^2/) < C. 

\Pi{t)\ Jn.{t) v ) Jn.it) 9 J 
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Let a and /3 be the two positive roots of the equation $( 2 ) = C. Then we have 


a < 


1 

rar 



~(s,y)dy, 

V 


1 

¥W\ 


[ Us,y)dy<(3. 

In i(t) 0 


These estimates imply (2.22). Finally we employ the mean value theorem to (2.22) to find cq(s, t ), bi(s , t) £ 
satisfying (2.23). The proof of the lemma is completed. □ 

We deduce a local representation of v in the next lemma by modifying Jiang’s argument in [6, 7] for 
fixed domains. To this end, we introduce the cutoff function ip z £ VF 1,00 (M) with parameter ^ £ R by 


<Pz(y) = < 


1, y<[z\+ 4, 

[z\ + 5 - y, [z] + 4 < y < [z] + 5, 

0, y > [ 2 ] + 5. 


(2.24) 


Lemma 2.5. Let (t,z) £ f It- Then we have 


j{t,y) = B.(t, y)A z (t) + — f V) ds 

d Jo B z (s,y)A z {s) 


B z (s,y)A z (s) 

for all t £ [0,t] and y £ I z (t) := (Y(t),oo) D ([ 2 ] — 1, [z] + 4), where 


B z {t,y) :=v 0 (y)expl- I (u 0 (£) - u(i,£)) ip z (£)d(, 
l M J y 

A z (t) := exp l-f f [Z]+5 - p) • 

[dJ0 J[z}+ 4 V V / J 


Proof. We multiply (2.19) 2 by ip z to get 


(<P*u) t = - p) - L>'z U- - p) ■ 

LVu / J v V v / 


(2.25) 

(2.26) 

(2.27) 

(2.28) 


Let (t,y) £ [0,r] x J z (r). Since y > Y(s) for each s £ [0,r], we have [0, r] x [y, 00 ) C f It- In view of the 
identity <p z (y) = 1 and (2.19)i, we integrate (2.28) over [0 ,t] x [y, 00 ) to get 


Pz(Q(u(t,£) - u 0 {£))d£ = - R [ Q[ S,V } ds 


r d(s,y) j [* [M +5 ( p _ i uy\ 

Jo v(s, y) S Jo J [z]+ 4 v M v ) ' 


This implies that for each t £ [0,r] 

1 


>{t,y) 


exp 


{-/ 

l M Jo 


vo(y) 

1 0{s,y) 


ds> = 


(s, y) J B z (t, y)A z (t)' 


(2.29) 


Multiplying (2.29) by and integrating the resulting identity over [0,t], we have 


exp 


l d Jo 


0{s,y) ^ 

v{s,y) 


ds )■ = 1 H-- 


Q{s,y) 


h Jo B z (s, y)A z {s) 


ds. 


We then plug this identity into (2.29) to obtain (2.25) and complete the proof of the lemma. □ 

The following lemma is devoted to showing the bounds of the specific volume v(t, z ) uniformly in the 
time t and the Lagrangian variable 2 . 

Lemma 2.6. If (2.6) holds for a sufficiently small eo > 0, then 


C 1 1 < v(t, z) < Ci for all (r, z) £ 0 T . 


(2.30) 
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Proof. Let (r, z) £ f It be arbitrary but fixed. The proof is divided into three steps. 

Step 1. In view of Cauchy’s inequality and (2.20), we have 

5 z (f,y)~l for all (t,y) £ [0, r] x I z (t). (2.31) 

Let 0 < s < t < t. Apply Cauchy’s inequality and Jensen’s inequality for the convex function 1/x 
(x > 0), (1-9), (2.20) and (2.22) to deduce 


t r[z} +5 


a 


< c 


< c 


L V 


p 


1+4 

t ^[z]+5 u 2 ^ i~[z ]+5 q ^[z]+ 5 rq 


s J[z}+ 4 

n [z]+5 

,t]+4 


vO 


v9 v6 


1+4 

R 


1+4 


[z]+5 Q 
:]+4 v 


(2.32) 


r* 

/■N+5 

/ inf 0(i',-) 

/ « 

Is (N+4.H+5) 

J[z]+4 


dt' 


<C + Ce 0 (t-s)-C~ 1 [ inf 9{t',-)dt', 

Js (W+4.W+5) 

For each 0 < t' < t, there exists y{t') £ [[z] + 4, [z\ + 5] such that 

= inf 6>(f',-), 

(W+4.N+5) 

Since [z] + 4 > [F(t')] + 2, we derive from the definition (2.21) that fl[ 2 ]+ 4 (f') = [[z] + 4, [z] + 5]. We then 
apply Holder’s and Cauchy’s inequalities to obtain from (1.9), Lemma 2.4 and (2.20) that 


a vd‘) a 

mat' 

M+4<y,t') " 


< 


S + fi[ z ] + 4 (i') 




ft fOO -= 

<mf 1 S{t-s)+ / -jf{t',C)df / v{t',f)df 

Js Jy(v ) vV Jn, r ,, A (t') 


dfdt' 


< C(t — s) + C 

< C(t — s) + C. 


Y(t') 
t poo 

Y(t') VV 


n W+ 4(t') 


dt' 


(2.33) 


Applying Jensen’s inequality to the convex function e x , we have from (2.23) and (2.33) that 


/ «, i n r f i ^ e i t '^) dt '= [ ex P (hi 9(t',y(t')))dt' 
J s (H+4.H+5) J s 

> (t — s) exp ^ J In 9{t ', y{t'))dt' 


> (t — s) exp 


t — s . 


rvO 


4(*',*') 


'-(t',£)d£,+ln9(t',b [z]+4 (t',t')) 


dt' 


> {t — s) exp — In Cq — 


t — s 


t i-yO 


s Jb [z j +4 (t',t') 


-f{t',Odfdt' 
/ t /\ w 


> 


f — s 


C 


■exp - 


C 

t—s 
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This implies 

for 0 < t — s < 1, 

—C _1 (t — s) for t — s > 1. 

Plugging (2.34) into (2.32) and taking eq > 0 small enough, we have for each s £ [0, t\ that 

ft f[z]+5 r 


/ inf 9(t',-)dt' < < 

J s ([*]+4,[*]+5) [-C - ( 


(2.34) 


// 

J S J \l 


M+4 


filly 


- P 


<C -C~\t- s). 


According to the definition (2.27), we then obtain 


0 < A z (t) < Ce ~ t/c , 


A *(t) 

A z {s) 


< Ce- {t ~ a)lc for all 0 < s < t < r. 


Step 2. Plugging (2.31) and (2.35) into (2.25), we infer 

f Tr\°( s ^y) ds ~ v (^y) < 1 + [ 0(s,y)e~^ds 
Jo A z{s) J o 


(2.35) 


(2.36) 


for all (t,y) £ [0,r] x J z (r). In light of the fundamental theorem of calculus, we infer from (1.9) and 
(2.22) that for y £ 7 z (r) and 0 < s < t < r, 


0{s,y)* ~ 9(s,b [z]+2 (s,T))z 
< [ 9-i\e y \(s,0d£ 

L 


< / e~2 

’I At) 

_ 1 

< 77ln 2 (5 + 


(s,tn+ / 9-?\d y \(s,Z)dl; 

I At) 


(2.37) 


f dl 

2 

f 



/ v9(s,£)d£ 

JiAP 


< m 2 2 <5 + sup u 2 (s, •) 

!z(t) 


r d 2 
11 AO Vd 


where we have used 6 [ z j + 2 (s, t) £ SA [^]+ 2 (t) C 7 z (t). Combine (2.37) with (2.23) and (2.6) to give 

d 2 

I 

A(r) 

and 


r d 2 

0(s, y)< 1 + sup v(s, •) / -4(s, £)d£ 

I At) J I z (t ) v9 


1 - C sup u(s, •) [ -^(s, £)d£ ^ 0(s, J/)- 

/*(->■) 7/,(r) v9 

We plug (2.38) into (2.36) to obtain 

/*£ r ^2 

v(t,y)<l + sup v(s,-) -^(s,£)d£ds. 

Jo i At) j I At) v 6 


(2.38) 

(2.39) 


Taking the supremum over / z (r) with respect to j/, we have 


P r d 2 

sup v(t,-) < 1 + / sup v{s,-) -^{s,Z)d£;ds. 

I At) Jo iAt) Jo. At) w 


(2.40) 


Applying Gronwall’s inequality to (2.40), we can deduce from (2.20) that 

sup v(t, •) < Ci for all t £ [ 0 , r], 


(2.41) 
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where C\ > 0 is some constant independent of t, r and z. Noting that z G J z (t), we deduce from (2.41) 
that v(t, z) < C i- Since (r, z) G f It is arbitrary, we conclude 

v(t, z) < Ci for all (r, z) G fir- (2.42) 


Step 3. On the other hand, in view of (2.22), (2.31) and (2.35), we integrate (2.25) on / z (r) with respect 
to y to find 


1 < 


Iz(t) 


v{t,y)dy < e 


-t/c 


+ 


Mi) 

A z(s) 


ds. 


Consequently, we have 




Jo M s ) 

Inserting (2.39), (2.42) and (2.43) into (2.36), we have 


ds> l-Ce~ t/c . 


v(t,y)> f ^lds-C ft Mt) 


/o 


A 2 (s) 


'0 


M s ) 


'!z(t) 


d 2 

^ ds 


> 1 - Ce~ t/c - C 

> 1 - C&-^ c - c 


Mi) 

Ms) 


Iz{j) 


d 2 

^ ds 


>l-Ce~ t/c -Ce~^-C / ~b d £,ds 


r d 2 f* r d 2 

/ ~^ d ^ ds - C / ~k d ^ S 

h m (r) v0 J | Jlz(r) Vd 2 

dl 


Iz(t) 


v6 2 


(2.43) 


(2.44) 


> 1 for all (t, y) G [T 0 , r] x J z (r), 

where To is a positive constant independent of t. In particular, the estimate (2.44) implies 

v(t, z ) > 1 f° r a ll T ^ T 0 , z > Y (r). (2.45) 

As in [12, 13], we can derive a positive lower bound for v, that is, 

v(t, z) > e~ Ct , for (t,z) G Ht- (2-46) 

Finally, we combine (2.46), (2.42) and (2.45) to get (2.30). This completes the proof. □ 

As a corollary of Lemma 2.6, we get the uniform bounds for the density p in the Eulerian coordinate. 
Corollary 2.7. If (2.6) holds for a sufficiently small eo > 0, then 

Cff 1 < p(t,x ) < Ci for all (t,x) G [0,T] x R+, (2.47) 

where the positive constant C\ depends solely on inf R+ (po,0o) and ||(<fo),V , o>$o)||i- 


2.4 Pointwise estimates of temperature 

In the following lemma, we employ the maximum principle to get the lower bound for the temperature, 
which does depend on the time t. 

Lemma 2.8. Suppose that (2.6) holds for a suitably small eo > 0. Then 


inf > 

R + 


inf R+ G(s, •) _ 

C 2 inf R+ 6(s,-)(t - s) + 1 


for 0 < s < t < T, 


(2.48) 


where the positive constant C 2 depends only upon inf R+ (po, Go) and ||(<(>o, ifo, i?o) II l • 
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Proof. It follows from (1.1 ) 3 that 9 satisfies 


Ot ~\~ U0 X @xx — 

CyP CyP 


2 P 

U x - U X 

ft 


> - 


P 2 _ R 2 p 
‘ fpc v p 4,pc v 


e 2 . 


Hence we deduce from (2.47) that 


6 t + u9 x - 9 XX + C^d 2 > 0. 

Cy P 


Let 0 := 9 - 0 with 0 := g 2 ■ We observe 


and 


©U=0,oo > 0, B|t= s > 0, 


0* + uO x - 0 X x + C 2 [9 + 0)0 — 9t + u9 x 9 XX + C20 2 > 0. 

CyP CyP 


Applying the weak maximum principle for the parabolic equation, we have that Q(t,x) > 0 for 0 < s < 


t <T and x G R+. This completes the proof of the lemma. 

Next we have the P 2 -norm in both time and space of d x . 

Lemma 2.9. If (2.6) holds for a sufficiently small positive constant eo, then 


□ 


sup / [ (j) 2 +d 2 +^]+ T f [(l + 0 + ^l + dl\<l. 

0 <i<TjR, Jo J Hj. 


r T 

/ 0 -/1K+ 


(2.49) 


Proof. We divide the proof into five steps. 

Step 1. First, for each t > 0 and a > 0, we denote 

fl' a {t) := {x G R+ : d(t,x) > a}. 

Then it follows from (2.8) and (2.47) that 

rt 1 fT 


sup 

0 <t<T 


f <t> 2 + f d 2 + / |<K S ,0)| 2 ds + f f [^l + d 2 x ] 

JR+ JR+\Q' a (t) JO Jo Jr+XCILU) 

r t 


< C(a) sup 

0<t<T 


J p£ + J 

Jr+ Jo 


) (s,0)ds 


0 J«. + \n’ a (t) 

T r , € 
0 JR + . 


■C(a) f f 
Jo J 1 j 


(2.50) 


X I X 

9 9 2 


< C(a). 


Step 2 . We now estimate the integral f 0 J n , ^ d 2 . For this purpose, we multiply (2.2 ) 3 by (d — 2)+ := 
max}!? — 2,0} and integrate the resulting identity over (0 ,t) x R + to obtain 

r t r rt 


% f P(d~ 2) 2 + + k[ [ d 2 +R f j p9M<>- 2)+ 

A J r + Jo Jn' 2 (s) Jo Jr+ 

= t / Po^o- 2) 2 + + f f h(d- 2) + + p[ t f 1%(d-2)+. 

z Jr+ Jo Jr+ Jo Jr+ 


(2.51) 


To estimate the last term of (2.51), we multiply (2.2 ) 2 by 2 ijj(d— 2 ) 4 . and integrate the resulting identity 
over (0, t) x R + to find 

[ if 2 p{d - 2 ) + + 2p ( f ifl(d- 2 )+ - [ ( pTp 2 (dt + ud x ) 

J r + Jo J r + Jo Jn' 2 (s) 


%Po(d 0 -2)+ + 2R / p9if x (d — 2)_|_ + 2R / / P 9f,d x 

R + Jo j R+ Jo Jn i(s) 

2 r( f ip(p9) x (d - 2) + - 2p f f ifip x d x + 2 f f gijj(d 
Jo J R+ Jo Jn’ 2 (s) Jo J R 


(2.52) 


■+ 
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Combining (2.52) and (2.51), we have from (2.2 ) 3 that 


L 2 

Cy 

.~2 


■p{9-2) 2 + + i/j 2 p{i}-2)+ 


0 J £l' 2 (s) 




^(*- 2 )+ 


0 ^R+ 


Po($o — 2)+ + ^opo^o — 2)+ + ^ J p , 

p=i 


where each term J v in the decomposition will be defined below. We now define and estimate 
terms in the decomposition. We first consider 


Ji ■= R f [ pOM#- 2)+ 

J 0 J R+ 


and 


In light of (2.47) and (2.8), we have 


Ji := 2 R 


pO^Vx- 


o Jn’ 2 (s) 


ip 2 + [ 0< f P £< 1. 

M_|_ 

From Cauchy’s inequality and (2.47), we obtain 

\Ji\ <*f f ^-2)+ + C(e) f f e\V-2) + 

JO J R+ Jo J R+ 

f [ i’l(v-2) + + c(e) f [ e&-i)l 

Jo Jr+ Jo Jr + 

V£(0-2) + + C(e) f sup(tf - 1) 

Jo R+ 


and 


< e 


\J*\< 


< e 


< e 


2 

+ > 


/O J Cl 2 (s) 

rt 


+ C(e) f f 
J0 Ji 2 ' 




1p 2 d 2 


j f #l + C(e)[ f V’ 2 ^-!) 

Jo Jq 2 (s) Jo J q' 2 (s) 

[ [ # 2 X + C(e) [ sup($ — 1) 

Jo Jn'(s) Jo R+ 


)%■ 


Let us define 


Jo ~ [ [ [h('V-2) + +c v 1 hilj 2 + 2gip('d-2) + \ . 

Jo JQ' 2 (s) 


It follows from (2.47) and the identity P — P = Rpd + R9(f> that 

Iffl ^ |wx|(</>,- 0 )|, \h\ < \{u x ,O x )\\((t>,ii,'d,i) x )\. 
Applying Cauchy’s inequality to J-$ yields 


\Js\< 


' 0 J ^2 (®) 
rt 


\{ux,O x )\\{<l),ip,i},ip x )\ [(tf- 2 )+(l + + ip 2 ] 


< 


0 JQ 2 (s) 


\{u x ,e x )\ [((/) 2 + ip 2 + V 2 ) + ipl + (d - 2) 2 + (l + ip 2 ) + ip 4 ] . 


(2.53) 

all the 


(2.54) 


(2.55) 


(2.56) 


(2.57) 


(2.58) 
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2 

+ > 


We obtain from (1.9) and (2.54) that 

[ [ |(uxA)|(tf - 2)^(1 + ip 2 ) < [ sup(tf-l) 

Jo J R + Jo R+ 

/ [ \(Ux,0xM 4 < [ W\\t°°\\{u x ,O x )\\ L i <s 

JO J R+ Jo 

Plugging (2.13), (2.59)-(2.60) into (2.58), we deduce from (2.50), (2.54) and Lemma 2.2 that 


) 2 + 


\Ja\ <$j [|^(s, 0)| 2 + \\{(/)x,ipx,^x){s)\\' 2 ] ds + J sup(i?-i) 

~ 1 + S IL +M| (7 + 7 )] + l ^“ 11 

<l+[ sup( , d - 1 ) 

Jo R+ 


4- 


Let us now consider the term 


which is trivially estimated by 


Ja ■■= 2 R 


pQ) x (fl - 2 )+, 


\M < / sup($ — 2 ) + ||' 0 ||||(p x ,# x )|| < / sup($ !)+■ 


to M+ 


to M + 


For the term 


Jo ■■= I f [flip 2 x c v V 2 - PtpxC v V 2 - 2 iri/jipx'$ x \ , 

Jo Jw 2 {s) 


we apply Cauchy’s inequality and (2.54) to deduce 


\Js\<ef ( € + C{e) 

Jo J Q' 2 (s) 


to J Q,' 2 (s) 

rt 


// 

Jo Jnu 


^V 2 + / / ^e 2 

O) 

t 


< 


< 


f f *l + C(e) [ f ^ 2 X + f [ if (0-1) 

Jo Jq' 2 (s) Jo Jq' 2 (s) Jo Jq 2 (s) 

ft € + C{e) f f + T sup($ - 1)7 

Jo Jnt(s) Jo JsiUs) Jo R+ 


0 J Cl' 2 (s) 

We finally consider 

Jo := / [ c~ x Kill'd xx . 

Jo J Cl 2 (s) 

In order to estimate f/e, we apply Lebesgue’s dominated convergence theorem to find 

r t 

I C,cn„/,2 

J XX 


Jo — lim [ [ ip v {f))ii 2 d x 

c v v-*o+ J q j R+ 

= — lim [ [ [~2tp v ('&)'ijjip x '& x - (pl('d)’ip 2 'dl 

C v I'->0 + J Q J R+ 

— lim [ [ 2 ip„('d)'ipip x 'd x 

C v V^O+ J q 


< 


< e 


f f € + C(e) j f 
Jo J R+ Jo J Rj- 


0 J R+ 


(2.59) 

(2.60) 


(2.61) 


(2.62) 


(2.63) 


(2.64) 
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where the approximate scheme y?„( 0 ) is defined by 


1 , 

(0-2)/v, 

0 , 


0 - 2 > v, 

Q <V-2 <v, 

V- 2 < 0 . 


<Pv(&) = < 

Plugging (2.55)-(2.56), (2.61)-(2.64) into (2.53), we get from (2.47) that 


< 1 + e 


Step 3. We obtain from (2.8) that 


/ (0-2)* + / / [0 *+^- 2) +] 

JR + J0 JO' 2 (s) 

f! € + C{e) f f + C'(c) f sup(^» - I)* . 

Jo lix Jo JRj- Jo r+ 


(2.65) 


+ / / [tfx + 2 ) + ] + / / 

«/0 Jfl'c,(s) J 0 ^R+\f 2 g(s) 


/0 «/f2g(s) 

</ 4 / [0® + '0x(0 — 2)+] + 1- 

./0 J Q' 2 (s) 


e 2 e 


( 2 . 66 ) 


Combining (2.66) and (2.65), and choosing e sufficiently small, we have 

f (0 - 2 ) 2 + + f f [€ + ^0] < 1 + C sup(tf - 1) 2 + + f f ^ l . 
Jrj. Jo Jro. Jo R+ Jo Jr+ 


(2.67) 


Step 4. To estimate the last term of (2.67), we multiply (2.2 ) 2 by ip 3 and then integrate the resulting 
identity over ( 0 , £) x R+ to have 

1 


pip A + 3p 


V’Vx - T / PoV’l 


0 


> 0 •/ M+ * «/ J&+ 

= 3 R I I ip 2 ip x 9(p + 3 R f f ip 2 ip x pi) + f f gtp 3 . 
I o Jr+ Jo Jr + Jo Jr+ 


( 2 . 68 ) 


From (2.50) and (2.54), we have 


/ f ip 2 ip x 9(p+ [ f ip 2 tp x p'd 

Jo JR+ Jo JR + \f 2 ^(s) 


(2.69) 


< 


L°°\\Vx 


ll^llL 2 (R + \ni,(s))] / IIV’xll" 

Jo 


We then apply Cauchy’s inequality to derive 


2„o2 


/ [ ip 2 ipxP# < e / / V’Vx + C'(e) / / 

Jo Jo' 2 (s) Jo Jr + Jo Jn 2 (s) 

ft r rt 

< e 

/o Jr + Jo r+ 


f [ iP 2 iP 2 x + C{e) f su p( 0 -l) 
Jo Jr, Jo R+ 


(2.70) 
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In view of (2.50), (2.54), (2.57) and (2.7), we have 


[ [ g^ 3 ^ / f \u x \ C0 4 + 0 4 ) 

JO J R+ Jo JR+ 

<sJ\\\mM 2 +m 2 m\ 2 ) 

<6 fu x r+5Mim^ f* [ M 

Jo Jo Jr + P 

<6 [ \\f x \\ 2 + 6M 2 mf 1 M 2 . 

Jo 

Plugging (2.69) -(2.71) into (2.68), and taking e sufficiently small, we derive from (2.6) that 


(2.71) 


/ 

J r+ Jo Jr 


rr x < i + 

0 J R+ Jo 


4>l+ [ sup(tf-l) 
Jo R+ 


2 

+ • 


(2.72) 


We note from (2.8) that 


< e f j (hp 2 + C( e ) [ f - 7 ? < e [ [ + C( e )- (2.73) 

Jo JRj. Jo Jr + a Jo Jr + 


Combination of (2.73) and (2.72) yields 


f </> 4 + f f (1 + V> 2 )^ < C(e) + e f f e^l+ f sup (0 - 1 ) 2 + . (2.74) 

Jro. Jo Jrj. Jo Jro. Jo r+ 


We plug (2.74) into (2.67) and choose e suitable small to find 


(0-2)1+^+ / [til+rt(l+e + il; 2 )\ <1 + / sup ( 0 - 1 ) 


/0 M + 


+ • 


(2.75) 


Step 5. It remains to estimate the last term of (2.75). According to the fundamental theorem of calculus, 
we have from (2.54) that 


[ SUP (0 - 1) + < [ [ \0 

Jo R+ Jo Jo; is) 


-i 2 


0 JQi(s) 


< 


< e 


< e 


tf 2 

U X 

e 


n'(a) ° J O'(a) 

f T f 0l+C(e)f T f 

Jo Jr+ Jo Jr h 

f T f € + C(e). 

Jo Jr + 


(2.76) 


€ 
e 2 


Plug (2.76) into (2.75) and choose e > 0 suitable small to obtain (2.49). This completes the proof of the 
lemma. □ 

We obtain the upper bound for the temperature uniformly in both time and space in the next lemma, 
by combining Lemma 2.9 and some desired uniform estimates on the spatial derivatives of 

Lemma 2.10. If (2.6) holds for a sufficiently small e > 0, then we have 

9(t,x) < Co for all (t,x) £ [0, T] x R + , 


sup ||( 0 ,^M)(*)lli + / 
0 <t<T Jo 


ll^(*)l| 2 + ll(^x,^)(t)||? 


dt < C\. 


(2.77) 

(2.78) 
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Proof. First, plugging (2.47) into (2.17), we deduce 


[ fa + [ [ 0<& < 1 + f [ 

J R. JO J R, Jo J Rj 


i9 2 

tb 2 + tf 2 + -2. 
rx w u x T 


< 1 , 


(2.79) 


where we employed (2.8) and (2.49) in the last inequality. 
Next, multiply (2.2 ) 2 by to derive 


2 




, 1 „ , 1.2 , ,>xx _ ( p - p )x„, 9,1. 

H - P - Txx Jpxx* 

2 p p p 


Integrating this last identity over (0,T) x R + , we obtain from (2.47) and Cauchy’s inequality that 


V’x-/ ^l(t,Q)dt + [ [ i/>. 

Jo Jo j R+ 


2 


< 1 + 


CP--P)x+S + ISxIV’x + IV>X 


^0 ^R+ 

Apply Sobolev’s inequality and (2.49) to obtain 


(2.80) 




l^l 3 < / lhMlhMI+ / lhMI*lliM* 


< 


[ T [ Tpl+C(e) [ T \\m 2 + lllM^ 
/o ./R+ J 0 

r T 


< 


C(e) 


1 + sup ll'i/’xll 3 
0 <t<T 


/ 

J Rj 




(2.81) 


•/0 ./R+ 

In view of the identity (P — P) x = R{0(j) x + + <j)6 x + ’dpx), (2.57), (2.79), (2.49) and (2.13), we derive 


' 0 J R+ 


< 


(P-P^+^ + l^ 


\(9(f> x ,'ijj x ,'dx)\ + \(Px, u x , 0 X )\ \(<p, ip, $)| 


Jo JR+ 

^ 1 + II^IU oo ([ 0 ,T]xR + )- 


We plug (2.81)-(2.82) into (2.80) to get 


(2.82) 


sup \\4>x(t)\\ 2 + / V’L ^ 1 + II 6 'IU-([o,t]xr+) + sup HV’xll 3 - 

o <t<T Jo JR + 0 <t<T 


Then Young’s inequality yields the estimate 


sup ||^x(i)||“+ f [ V’xx ~ 1 + l|0||i“([O,r]xR + )- 
0 <t<T Jo j R-L 


(2.83) 


Next, multiply (2.2 ) 3 by and integrate the resulting identity over (0,T) x R + to have 


f f= c ij c + /7 

J R+ Jo J R+ P z JR+ JO JRj. 


/ 2 

c v wd x — p— + R9ip x —- 
P P J 


u xx-> 
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which combined with (2.47) implies 

[ €+ f [ e<i + /7 + 11^11 l~€ + o 2 € + h2 ) 

J r+ Jo Jr-l. Jo Jr+ 


10 j IK.|- 

<1+ / T (i + ||V , llll'0x||)||^x|| 2 + / T ||^ii 3 ii^ii 
Jo Jo 


r T r r T 

1. 2 , / 


II^IU“([o,r]xR+) / / 0$ 

0 t'M-)- J 0 


/i 2 . 


By (2.57), (2.13) and (2.83), we have f Q / R+ /i 2 < 1, and 

/ II0*|| 3 ||<MI < sup ||Vg| 2 f (ll^|| 2 +ll^xxll 2 ) < 1 +||6'|||=o ([0T]xR )■ 
Jo 0 <t<T Jo 

In light of (2.79), (2.83) and (2.49), we then obtain 


[ $1+ [ [ ~ 1 + l|0|l!°°([O,T]xR + )' 

J R_i_ J 0 J R-i_ 


*+ JO JR+ 

Finally, it follows from (2.49) and (2.85) that 

11^ — ^lll“([0,T]xR + ) ~ SU P 11^(01111^(^)11 < 1 + ||0||l”([O,T]xR+)- 
0 <t<T 

This implies (2.77) by virtue of Cauchy’s inequality. Combine (2.80), (2.83) and (2.85) to give 


sup [ [<(> 2 + 0 2 + tf 2 ] + / / + + 

0<t<T Jm., Jo Jr, 


0 JR+ 


which together with (2.49) yields (2.78). This completes the proof of the lemma. 


(2.84) 


(2.85) 


□ 


3 Proof of Theorem 1 

In this section, we prove Theorem 1 in six steps by employing the continuation argument. 

Step 1. Set T\ = 128A3 4 C|. We choose positive constants II, and A,; ( i = 1,2,3) such that 
||(<£o, 0o,$o)||i < n and 

Ai < po(x) < Ai, A2 < 9 0 (x) < A 2 , A3 < Q(x) < A 3 for all x £ M+. 

Applying Proposition 2.1, we see that there exists a constant 0 < C < min{Ti, T 0 (Ai, A 2 , II)} such that 

5A1,2Ai; f. 


the problem (2.2) has a unique solution (</>, 0,t?) £ X(0,fi; ^Ai,2Ai; ^A 2 ,2A 2 ). 


Let 0 < 5 < <5i with 


" 1 -y,2Ai, -y,2A 2 ) = e 0 . 


Then we can apply Lemmas 2.6, 2.8 and 2.10 with T = ii to obtain that the local solution (</>, 0,$) 
constructed above satisfies for each t £ [0,ti] that 


6(t,x) > 


A 2 


C 2 A 2 Ti + 1 


for all x £ R+, 


and 


C x 1 < p(t,x ) < Ci, 9(t, x) < Co for all x £ M+, 

||(0,0,i?)(t)|| 2 + f \\\VeMs )\\ 2 + \\(^Vz)(s)\\l 

Jo L 


ds < C\. 


(3.1) 


(3.2) 
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Step 2. If we take (0, ip, t ?)(£1 , •) as the initial data, we can apply Proposition 2.1 and extend the local 
solution (0, ip, &) to the time interval [0, t\ + £ 2 ] with £2 < min{Ti — £ 1 , g 2 a^+i 1 Cl)}- Moreover, 

we have 

2 k - p{t • x) - 2Cu w,u) s e(i ' 11 s 2Cs 

for all (£, x) £ [£i,£i + £ 2 ] x R+- Take 0 < 5 < min{<Ji, £> 2 } with 

2Sr' 2Cl -2(C 2 A2T 1 + l)’ 2C3 ) &=, “- 

Then we can employ Lemmas 2.6, 2.8 and 2.10 with T = t\ +£2 that the local solution (0, 0,i?) satisfies 
(3.1) and (3.2) for each £ £ [0 ,£1 +£ 2 ]. 

Step 3. We repeat the argument in Step 2, to extend our solution (0, ip, d) to the time interval [0, £1 +£2 + 
£ 3 ] with £3 < min{Ti — (£1 + £ 2 ), T 0 (-^, C 4 )}. Assume that 0 < 5 < min{<5i, £ 2 }. Continuing, 

after finitely many steps we construct the unique solution (0,0,$) existing on [0,Xi] and satisfying (3.1) 
and (3.2) for each £ £ [0,Ti]. 

Step 4. Since T\ > 128A^" 4 C4 and 


sup 

0<t<Ti 


pTi 


ITi 


WMmidt < cl 


we can find a t' 0 £ [Ti/2,Ti] such that 


X 2 

A 3 


Sobolev’s inequality yields 


llf«)ll<C4, K(OII<^- 


IWto)lk»»<^(£ , o)||’||t?x(£( ) )|| i <-^, 


and so 


0 (t' o ,x) > 0(x) - ||$(£o)||l“ > y for all x £ K+. 


We note here that 

||(0,0,$)(^o)|| 1 < C 4 , Cf 1 < p{t' 0 ,x) < Ci, 0(t' o ,x)<C 3 for all x £ R+. 

Now we apply Proposition 2.1 again by taking (0, if), i ?)(£q, •) as the initial data. Then we see that the 
solution (0, ip, 1 }) exists on [t' 0 ,t ' 0 + £j] with t[ < min{Ti, To(^-, ^ 3 , C 4 )} and satisfies 

2^ < $(C *) < 2Ci, ^ < 0(£, x) < 2C 3 
for all (£, x) £ [£q, £q + £ 4 ] x M + . If we take 0 < 8 < min{<5i, 82 , <$ 3 } with 

A.,2C 1 ,^,2C 3 )ia=e 0 , 

then we can deduce from Lemmas 2.6, 2.8 and 2.10 with T = t' Q + 1 [ that for each time £ £ [£q,£q + £)], 
the local solution (0,0,$) satisfies (3.2) and 


0 {t,x) > 


inf R+ 0 (t' o , •) 


> 


A 3 


for all x G M+. 


(3.3) 


C 2 inf R+ £»(£(,, -)Ti + 1 “ C 2 A 3 Ti + 2 

Step 5. Next if we take (0, ip, $)(£(, + t[, •) as the initial data, we apply Proposition 2.1 and construct the 
solution (0,0,i?) existing on the time interval [0, t' 0 +t[ +t' 2 ) with t' 2 < min{Ti —1 [, Tb(^v-, c . 2A ^ i+2 , C 4 )} 
and satisfying 

2 k S '’ <M) S 2C ” 2(CiA^T2)^ <! '") S2C3 
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for all (t, x ) £ [t ' 0 + t[,t ' 0 + t[ + t' 2 \ x K + . Let 0 < S < min{( 5 i, S 2 , ^3,^4} with 

“ ( 2 C[ l2Cl ' 2(C 2 X 3 Ti + 2 )’ 2 ° 3 ) Sa = £ °' 

Then we infer from Lemmas 2.6, 2.8 and 2.10 with T = tg + t\ + 1 / 2 that the local solution ( 4 i, if, d) satisfies 
(3.3) and (3.2) for each t £ [t' 0 , t ' 0 + t\ + t' 2 \. By assuming 0 < <5 < min{<5i, <J 2 , <^3, ^4}, we can repeatedly 
apply the argument above to extend the local solution to the time interval [0 ,£q + Tf]. Furthermore, we 
deduce that (3.3) and (3.2) hold for each t £ [t' 0 , t' 0 + T±\. In view of t' 0 + Ti > |Ti, we have shown that 
the problem (2.2) admits a unique solution ((f>,ip,d) on [0, | Tf\. 

Step 6. We take 0 < <5 < min{i5i, S 2 , S3, ^4}. As in Steps 4 and 5, we can find tg £ [t' 0 + T±/ 2 , t' 0 + Ti] such 
that the problem (2.2) admits a unique solution (cf>, d) on [0,tg +Ti], which satisfies (3.3) and (3.2) 
for each t £ [t' 0 ,tg + Ti]. Since tg +T\ > t' 0 + | T\ > 2Ti, we have extended the local solution (</>,'0,'d) 
to [0, 2Ti]. Repeating the above procedure, we can then extend the solution (<^>, “0, i?) step by step to a 
global one provided that 5 < min{<$i, S 2 , S3, ^4}. Choosing c\ = min{<5i, 62 , S3, A4}, we then derive that 
the problem (2.2) has a unique solution satisfying (3.2) and 


inf 6 (t , •) > min 
R+ 


A2 A3 i 

C 2 X 2 Ti + 1 C2A3T1 + 2 J 


for each t £ [0,00). 

Therefore, we can find constant C5 depending only on infR + (po,#o) and ||(</>Oj'0O)'$o)||i such that 

pOO 

sup \\(<l>,ip,'&)(t)\\l+ [||<Mi)|| 2 + ||(V>*,tf*)(£)||?] dt< Cl, 

0<t<oo J 0 

from which the large-time behavior (1.13) follows in a standard argument. This completes the proof of 
Theorem 1. 
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